Electric field suppression of ultracold confined chemical rates 
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We consider uhracold colhsions of polar molecules confined in a one dimensional optical lattice. 
Using a quantum scattering formalism and a frame transformation method, we calculate elastic and 
chemical quenching rate constants for fermionic molecules. Taking '"'K^^Rb molecules as a proto- 
type, we find that the rate of quenching collisions is enhanced at zero electric field as the confinement 
is increased, but that this rate is suppressed when the electric field is turned on. For molecules with 
500 nK of collision energy, for realistic molecular densities, and for achievable experimental electric 
fields and trap confinements, we predict lifetimes of KRb molecules of 1 s. We find a ratio of elastic 
to quenching collision rates of about 100, which may be sufficient to achieve efficient experimental 
evaporative cooling of polar KRb molecules. 
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The recent achievement of ultracold polar 
molecules (J-Q has opened tremendous perspec- 
tives to the ultracold community. The quantum states 
of the molecules can now be addressed experimen- 
tally 0, 3. Beyond the idea of creating molecular 
Bose-Einstein condensates or degenerate Fermi gases 
of polar molecules, these "quantum-state controlled" 
molecules have applications in ultracold chemistry [H, , 
condensed matter and many-body physics 0, Q , quan- 
tum information [T^j or in precision measurement 
and test of fundamental laws These applications 

crucially depend on the collisional properties of the 
polar molecules in the presence of an electric field. To 
get a stable gas of molecules, elastic processes should 
dominate over quenching (including inelastic and/or 
reactive) processes. This will be the case for certain 
molecules for which chemical reactions are energetically 
forbidden in their absolute ground state. For others, 
as KRb molecules, chemical reactions such as KRb -I- 
KRb — >■ K2 -l- Rb2 are energe tically allowed 
and evaporative cooling might be difficult to achieve. 
Chemical rates increase as the sixth power of the dipole 
moment induced by an electric field, due to head-to-tail 
collisions of the polar molecules [ij, [l^ . On the other 
hand, these collisions may be suppressed in a one 
dimensional (ID) optical lattice @, [l6| that confines the 
molecules in planes perpendicular to the axis of their 
polarization, so that the dipoles mutually repel each 
other. 

In this Letter, we extend the quantum formalism 
used in Ref. to describe ultracold collisions of po- 
lar molecules in an electric field confined in a ID optical 
lattice. We apply the theory to predict the elastic and 
chemical rate of confined KRb -f- KRb K2 + W02 reac- 
tions for which particular attention is devoted 0, 0, H, [3 
and for which theoretical predictions are needed. We 
show that the electric field suppression of confined chem- 
ical rates can help to achieve efficient evaporative cooling 
of such molecules. We assume that the wells of the ID 
optical lattice are deep enough so that the molecules can- 
not tunnel from site to site and can be approximated by 
independent harmonic oscillator traps. In the following, 
quantities are expressed in S.I. units, unless explicitly 



stated otherwise. Atomic units (a.u.) are obtained by 
setting h = Atteq = 1. 

We consider pairs of bosonic or fermionic polar 
molecules of mass m and reduced mass fi = m/2 con- 
fined in a harmonic oscillator trapping potential Vho = 
jiuP' which allows the molecules to collide in a two 
dimensional (2D) configuration space. The electric field 
is applied along the z direction, perpendicular to the 
planes of free motion in the trap. The Hamiltonian in 
the relative coordinate R is given by 



H = V" 
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ho 



dW 



dd 
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The long-range interactions are represented by the van 
dcr Waals potential Kdw = —Cq/R^ and the dipole- 
dipolc interaction Vdd — d^(l — 3cos^ 0)/(47reo R'^), where 
d represents the effective dipole moment of the polar 
molecule induced by the electric field Moreover, we 
represent chemical reactions via the non-Hermitian, ab- 
sorbing potential Kibs which has been chosen to success- 
fully describe three dimensional (3D) chemical rates of 
KRb molecules in an electric field, measured in Ref. [l^ . 

Scattering wave functions 4* are most naturally de- 
scribed at large intermolecular spacing in cylindrical co- 
ordinates R ~ (p, z, if) in accordance with the sym- 
metry of the confined trapping potential. The interac- 
tions between molecules when they are closer together, 
however, are better described in spherical coordinates 
R = {R,6,(p). In our numerical calculations, we will 
therefore use the appropriate coordinate system and then 
weld the wave functions together at a suitable matching 
distance. Note that the azimuthal angle (p is common to 
both coordinate systems, and in fact the Hamiltonian is 
independent of ip. Therefore, the quantum number M, 
representing the azimuthal projection of the orbital an- 
gular momentum, is rigorously conserved. Note also that 
in the limit of strong confinement, 6 is restricted to values 
near 6 « 7r/2, in which case the dipolar interaction is re- 
pulsive. This is the primary principle behind the electric 
suppression of collisions. 

When the molecules are close together, we solve the 
coupled-channel equations of motion in spherical coordi- 
nates using the diabatic- by-sector method [13, [13 • We 
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divide the complete range of R into sectors labeled by 
an index p. In each sector, we expand the M-dcpcndcnt 
wave function of an initial channel i as 



1 

R 



Y^xf{Rp;0)Ffl{R,-R). 



(2) 



The adiabatic functions xf {Rp'jf^) arc those that diago- 
nalize the Af— dependent angular part H^\R,d) of the 
Hamiltonian in Eq. ([1]) at the fixed radius R = Rp, mid- 
dle of the sector p. The associated adiabatic energy is 
ej(i?p), which converges to a harmonic oscillator energy 
En of a state n of the trap at large Rp. The radial func- 
tions Fjf{Rp\R), where / — 1,2,... represents an arbi- 
trary final channel of the system, are determined within 
each sector according to the diabatic equations of motion 



2fi dR^ 



^1 FfiiRp'^R) 

+ U}^ {Rp ; R) Ffl {Rp ; i?) = (3) 



where 



UfARp] R) 



xf{Rp;e)n^\R,e)xf{Rp; 



ain0d0. (4) 



E is the total energy of the system. To solve Eq. ([3]) we 
employ the method of the propagation of the logarithmic 
derivative matrix = {F^y^F^')' of Johnson ^ 
up to a suitable matching radius Rm- 

If Rjn is sufficiently large, then the only potential en- 
ergy of any significance is the trap confinement potential 
Viio. At this point the wave function is more conveniently 
described in cylindrical coordinates. Therefore for an ini- 
tial state Tii of the trap, we expand the Af -dependent 
wave function as 



M 



(5) 



where (z) are normalized harmonic oscillator func- 
tions in z corresponding to a state n/ of the trap, and 
represent the functions p^/^ xY / R in the asymptotic 
limit. The radial functions G^^ serve to define the reac- 
tance matrix K^^ via G*^.„^(p) oc p^/^ JM{kn;P) SnfUi + 

Jm,Nm are Bessel functions, 
the modified Bessel functions has 
f = y'2 //(£■- £nf)/h repre- 
sents the wave-vector of the relative motion of the state 
Uf. The K^'^ matrix is found by a matching procedure 
from the spherical wave function that captures the short- 
range physics and the cylindrical wave function that cap- 
tures the asymptotic boundary conditions. This is done 
by equating Eq. to Eq. (0) at a constant radius 
R = Rm [iTl . [isj , and taking into account the one-to-one 



For closed channels 
to be used instead 



interaction 


CtvdW.dd.ho 


(ao) i5vdw,dd,ho (/iK) 


Ce = 21000 ao 


264 


19.6 


d = 0.1 D 


179 


85.2 


a = U.o JJ 




i.UD 


d = 0.566 D 


5734 


0.083 


= 50 kHz 


1069 


2.4 


V = 1000 kHz 


239 


48 



TABLE I: Characteristic lengths and energies of the different 
interactions involved in the KRb 4- KRb collision. 1 D = 
1 Debye = 3.336 10"^° C m. 



correspondence between the short-range adiabatic states 
I, / and their long-range counterparts ni,nf. 

The K^^ matrix in turn determines the scattering ma- 
trix = {I + iK^^)^^{I - iK'^'), where / represents a 
diagonal unit matrix. The cross sections for elastic and 
quenching collisions are given by [20l - l22j 



h 
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niUi I 



X A 



(6) 



(7) 



X V, 



and the rate coeffici ent is given by /C^'^?" 
where v — y^2 Ec/fJ. is the collision velocity. The factor 
A represents symmetrization requirements for identical 
particles: A = 1,2 according as the particles are distin- 
guishable or indistinguishable. To compare with exper- 
imental results, one should average the rate coefficients 
over a Maxwell-Boltzmann distribution of the velocity v, 
but we do not perform this average here. Finally, due 
to exchange symmetry of indistinguishable molecules, 
the quantum numbers in the asymptotic representation 
Eq. ([5]) satisfy the relation 



(-1)^ 



-M 



= 7 



(8) 
-1 for 



with 7 = -|-1 for bosonic molecules and 7 
fcrmionic molecules. 

For concreteness, we now consider collisions of 
fermionic ''^K^^Rb molecules, prepared in indistinguish- 
able internal states. We use the actual mass and per- 
manent dipole moment of these molecules, and set the 
coefficient of their van der Waals interaction to Cg = 
21000 a.u. [lil (1 a.u. =1 E^a^ where Eh is the Hartree 
energy and ao is the Bohr radius). We consider only 
molecules initially confined to the ground state of the 
harmonic oscillator, with rii = Uf ~ 0, and we take the 
collision energy to be Ec = 500 uK, relevant in the on- 
going KRb experiment. Under these circumstances, we 
find converged results at a matching radius Rm = 10^^ 
flQ. As the pair of fermionic molecules is in the ground 
harmonic oscillator state n = 0, only odd values of M are 
allowed in Eq. ((S]). This circumstance removes the occur- 
rence of undesired head-to-tail collisions of the molecules 
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in an electric field at ultralow energy. At Ec = 500 nK, 
we find that scattering is largely dominated by the single 
partial wave M ~ I, similar to the 3D case [3,[il. De- 
spite the fact that the molecules are in the ground state 
of the trap, we need to employ one asymptotically open 
channel (n = 0) and three asymptotically closed channels 
(n = 2, 4, 6) to converge the results to 10 %. 

The characteristic lengths (a) and energies (E) of the 
different interactions involved in the chemical process 
are presented in Table [J They are given for the vari- 
ous interactions by a^dw = (2 /J. Cg/^i^)^/'* and E^^w = 
H"^/ (2 /i a^dw) '^^^ Waals interaction, add(d) = 

^.(P/h^ and Edd{d) — h^/{y?'d'^) for the dipole-dipole 
interaction, and a\io{v) = ^ h/{2TTv fi) = h/{uj fi) and 
Eho{j^) = fi2'Ki' = huj for the harmonic oscillator trap. 
This is useful for characterizing the different regimes in- 
volved in the chemical process. When Ovdw or Odd < aho, 
intermolecular forces take place where the confinement is 
small and the collisions are effectively 3D. When a^o < 
flvdw or add, the reverse is true and the collisions are ef- 
fectively 2D. For fermions in n = and when an electric 
field is applied, the molecules meet primarily side-by-side 
and repel each other. For weakly polarized molecules, 
the 3D and 2D limits are realized for trap frequencies 
of 50 kHz and 1000 kHz, respectively (see Table For 
stronger polarized molecules at d > 0.3 D, the 2D limit 
is reached for both confinements as a^o < add- We con- 
sider the two confinements in the following. We note 
that the present theoretical formalism can treat both 3D 
and 2D limits in an electric field. This is in contrast 
with former theoretical studies [l^, [l^] which do not 
describe the electric field dependence and the 2D limit 
where Cho < Ovdw or add- 

We present in Fig. [1] the elastic and quenching rate 
coefficient as a function of the induced dipole moment 
d, for a trap frequency of i/ = 50 kHz (upper panel) 
and ly = 1000 kHz (lower panel). The elastic rates 
(blue lines) have the same trend for the two different 
confinements, increasing with the dipole moment. For 
d > 0.2 D, the elastic rates converge to a semi-classical 
formulation (thin blue lines) of pure 2D dipolar scatter- 
ing [2^, which increases linearly with d. For d < 0.2 D, 
the elastic rates are in better agreement with a pure 
2D Born approximation (dashed blue lines) that scales 
as d"* [25|. For d ~ D, the elastic rates take a fi- 
nite value which is determined by an unknown scattering 
phase shift which depends on the short-range potential 
of KRb-KRb. We note that the overall elastic rate is in 
better agreement with the pure 2D estimations for the 
larger confinement ly = 1000 kHz than for the smaller 
confinement 1/ = 50 kHz. as one expects. 

In contrast, the quenching rates (red lines) have a dif- 
ferent trend with dipole moment for the two different con- 
finement strengths. For = 50 kHz (upper panel), the 
quenching rate first decreases and then increases again 
as a function of d. This behavior already contrasts with 
the quenching rate in 3D collisions, which increases as 
d^ in the absence of z-confinement [iJ: [3] - The crucial 
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FIG. 1: (Color online) Collision of '*°K*''Rb molecules in the 
ground state of a confining trap. The elastic (blue) and 
quenching (red) rate coefficients are plotted as a function 
of the dipole moment d for f = 50 kHz (upper panel) and 
v = 1000 kHz (lower panel), for a fixed collision energy of 
Ec — 500 nK. A pure 2D semi-classical (SC) formula (thin 
blue line) and a pure 2D Born approximation (BA) formula 
(dashed blue line) has also been reported from Ticknor [25|] . 



difference comes from the fact that destructive head-to- 
tail collisions are removed by the confined geometry for 
fermions in the ground state n = (see Eq. ([8|)). Un- 
der even greater confinement (lower panel) the quench- 
ing rate continues to decrease with increasing dipole mo- 
ment [1^ nil , illustrating the electric field suppression of 
confined chemical rates. 

To better understand the qualitative trend of the 
quenching rates, we present in Fig. [2] the height 14 of 
the effective potential energy barrier corresponding to the 
lowest adiabatic energy curve ef^i{Rp), as a function of 
d for 1/ = 0, 50, 1000 kHz. We follow the qualitative argu- 
ments given in Ref. [isj that the behavior of the quench- 
ing rates are suppressed by the need for the molecules 
to tunnel through this barrier to the region of chemi- 
cal reactivity. For 1/ = 50 kHz, the barrier increases for 
small dipole moments, since the dipole interaction is re- 
pulsive for M = 1. However, anisotropy of the interaction 
couples different channels together. Therefore, at higher 
dipole moments, the lowest adiabatic curve is repelled 
more strongly from the others, ultimately lowering the 
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FIG. 2: (Color online) Barrier heights Vt as a function of the 
dipole moment d for u = kHz (black), v — 50 kHz (red) and 
ly = 1000 kHz (green). 



barrier height and increasing the quenching rate. For the 
stronger confinement v = 1000 kHz, the barrier height 
continues to grow as the dipole moment increases, lead- 
ing to continued suppression of quenching, by four orders 
of magnitude over the range of dipole moment shown. 
There are two ways to see this ongoing suppression of 
quenching. First, since aho(1000 kHz) < Ovdw or add, 
the molecules exert strong dipolar forces on each other at 
long range where they are still described by harmonic os- 
cillator states. Thus the interaction remains overwhelm- 
ingly repulsive and the molecules do not get close enough 
to react. Alternatively, we note that the trap confine- 
ment is tight enough so that the spacing between the 
harmonic trap levels is large, and therefore the repulsion 
between adiabatic channels is small. The incident adi- 
abatic channel is therefore less likely to be lowered due 



to the presence of other locally open adiabatic channels 
at short range. As a result, the barrier height continues 
to increase as d increases, driven solely by the increasing 
dipolar repulsion. 

The calculations show another interesting trend as con- 
finement is increased. Namely, at zero dipole moment 
the quenching rate is higher for the tighter confinement 
1/ = 1000 kHz. This is consistent with the barrier height 
Vb being smaller for the higher trap frequency (Fig. [2]). 
This difference arises from the fact that the collision en- 
ergy is measured relative to the asymptotic energy of the 
potential, which is the zero point energy of the confin- 
ing potential. The tighter trap has a higher zero-point 
energy (Tab.[l|, hence the apparent barrier is lower. 

Finally, Fig.[T]shows that for an experimentally achiev- 
able frequency trap of = 50 kHz (upper panel), we 
predict a loss rate of ^ 10^^ cm^ s~^ per molecule for 
the maximum dipole moment d = 0.2 D achieved so far 
in the KRb experiment. For a realistic planar density 
of molecules ^ 10^ cm~^, this corresponds to molecu- 
lar lifetimes of ^ 1 s which is 100 times longer than if 
the molecules were not confined. More important is the 
number of elastic collisions per chemical reaction, given 
by the ratio of elastic to quenching rates. At d ~ 0.2 D 
for = 50 kHz, we predict ^ 100 elastic collisions per 
chemical reaction. This ratio may be sufficiently high to 
achieve efficient evaporative cooling of KRb molecules in 
a ID optical lattice. 
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